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The three band p — d model of strongly correlated electrons interacting with optical phonon 
via diagonal and off-diagonal electron-phonon interaction is considered within cluster perturbation 
theory. At first step the exact diagonalization of the Hamiltonian of Cu 04 cluster results in the 
construction of local polaronic eigenstates |p) with hole numbers = 0,1, 2 per unit cell. The inter 
cluster hoppings and interactions are exactly written in terms of Hubbard operators = |p) {q\ 
determined within the multielectron polaronic eigenstates \p) at site f. The Fermi type single elec¬ 
tron quasiparticle dispersion and spectral weight are calculated for the undoped antiferromagnetic 
parent insulator like La 2 Cu 04 . The quasiparticle dispersion of Hubbard polarons is determined by 
a hybridization of the several Hubbard subbands with local Franck-Condon resonances. For small 
electron-phonon interaction the conductivity band is stronger renormalized then the valence band. 
Nevertheless for large electron-phonon interaction both bands are strongly renormalized with quasi¬ 
particle localization. Effect of partial compensation of diagonal and off-diagonal electron-phonon 
interaction at intermediate coupling is found. 

PACS numbers: 71.38.-k, 74.72.-h, 71.10.Hf 

I. INTRODUCTION 

The electronic structure of the undoped cuprates which 
are parent compounds for the high temperature super¬ 
conductors is one of the most important problem for 
understanding the physics of the high-Tc superconduc- 
tivityi. Strong electron correlations (SEC) are known 
to result in the Mott insulating state. The polaron for¬ 
mation and strong electron-phonon interaction (EPI) are 
expected from the ionic nature of the undoped cuprates. 

The interplay between EPI and SEC is a key to resolve 
the quantum dynamics of the doped holes or electrons 
into the high-Tc cuprates. 

Experimental indications for strong EPI have been 
provided by angle-resolved photoemission spectroscopy 
(ARPES) studies of the undoped Ca2Cu02Cl2 where 
broad Gaussian spectral features are found and inter¬ 
preted in terms of Franck-Condon processes and polaron 
physicsi^i^ The large isotope effect on Tc for underdoped 
cuprates and on the superfluid density at the optimal 
doping also suggests the importance of the EPl4i^ In ad¬ 
dition the ARPES data in doped cuprates have revealed 
the electron dispersion strongly renormalized by EPIji 
The interplay between charge, spin and lattice degrees 
of freedom has been studied theoretically by different 
approaches: the diagrammatic Monte Carlo method^ 
the exact diagonalization within the limited functional 
spaced the dynamical mean-field calculations within the 
Hubbard-Holstein model^^^ii and the determinant quan¬ 
tum Monte Carlo method/^ see the review paper 

To consider both the local effects of the strong EPI 
as a set of the Franck-Condon resonancesi^ and the 
electron dispersion in the infinite lattice we present in 
this paper the polaronic version of the multielectron 
generalized tight-binding (GTB) approach to the quasi¬ 
particle band structure in SEC materials. The gener¬ 


alized tight-binding methodic and its ab-initio version 
LDA-l-GTBi^ have been developed previously to study 
the electronic structure of the cuprates, cobaltites and 
manganites within the multiband Hubbard model. The 
polaronic version of the GTB method (p-GTB) as well 
as the initial GTB is an example of cluster perturbation 
theory] The general idea of the GTB is the follow¬ 
ing: the initial multiband Hubbard-like Hamiltonian is 
written as the sum of the intracell part and the in¬ 
tercell hoppings and interactions Hcc- The exact diag¬ 
onalization of the He provides a set of local multielec¬ 
tron eigenstates {|p)} and the intracell Hubbard opera¬ 
tors = Ip) (g|i^ Hence for, all two-site intercluster 
contributions in Hcc will be given by a bi-linear products 
of Hubbard operators of different unit cells. In other 
words, any multiband Hubbard like model in GTB looks 
like the original Hubbard model. The only difference is 
the size of the matrix |p) (g|. In the original Hubbard 
model with 4 local states it is the 4x4 matrix. In the 
5-band p — d modeli^ with number of local electronic 
states Ne Ri 100 the \p) (gj matrix size is Ng x Ng. In 
polaronic approacb^d^ with number of coupled phonons 
Nph > 10 the size of \p) {p\ matrix may be 1000 x 1000. 
Nevertheless for the low energy physics only a small part 
of multielectron states is relevant that is very important 
for the p-GTB approach below. 

Besides the local states of Cu d-holes and O p- 
holes and local Coulomb interactions Ud, Up, and Vpd 
we include the electron-phonon interaction with optical 
phonon mode. Exact diagonalization of the Hamiltonian 
Hg results in the local polaronic states {|p)} defined for 
different number of holes per unit cell nn = 0,1,2. For 
the three band p — d model2£i2i with EPI the same pro¬ 
cedure has been carried out in the Ref. (2^. Within the 
GTB perturbation treatment of the intercluster Hamil¬ 
tonian Hgg we have calculated the dispersion of Hub- 
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bard polarons that is formed by hybridization of the 
Hubbard subbands with the local Franck-Condon res¬ 
onances. To distinguish standard polaron for non- or 
weakly correlated electrons from polaron in strongly cor¬ 
related electronic systems we introduce the notion of 
Hubbard polarons for the latter case. The first doped 
hole or electron are delocalized at the weak EPI and be¬ 
comes localized at the strong EPI. Effects of the diago¬ 
nal (Hubbard-Holstein model) and off-diagonal EPI have 
been discussed. Partial compensation of these two EPI 
has been found at the intermediate coupling strength. 

The organization of this work is as follows. In the next 
section we present the three band p — d model Hamilto¬ 
nian with additional bare phonons and EPI and separate 
the total Hamiltonian to a sum of individual unit cell part 
He and the intercell contribution Hec- Section nni clar¬ 
ifies briefly the general ideas of polaronic version of the 
GTB formalism. In section llVl we discuss the structure of 
local Hilbert space with multiphonon and multielectron 
eigenstates with relevant for cuprates quantum numbers 
of holes per cell Uh = 0,1, 2. Sections IVl and IVll contain 


respectively the calculated polaronic band structure and 
spectral function for undoped cuprates. Discussion of the 
results is presented in Section IVHl 

II. THE MINIMAL MODEL 

Among the different groups of high-temperature super¬ 
conductors the copper oxides can be classified as the most 
correlated ones. To consider interplay between strong 
EPI and strong electron correlations in these compounds 
we start from the minimal semirealistic model which de¬ 
scribes holes in the CuO-plane interacting with the lon¬ 
gitudinal optical vibrational mode: 

H = Hpd + Hph + Hepi- (1) 

The Hamiltonian Hpd of three-band p — d model^Si^i in¬ 
cludes only orbitals that are most essential for the low- 
energy physics of cuprates, i.e. bonding Cn-d^'^-y'^ and 
0-Px,Py ones: 


Hpd 


(Erf - ^^) +Yi^p- d) P(a)gaP(a)ga 

f,cr Q;,g,cT 

H ip^ d {4aPic)sa + ^-C-) + ^ I] ^ 

Oijf^g.O' f,{7 


Here indexes f and g run through positions of copper 
and oxygen plane atomic orbitals so that g = f -f and 
index I enumerates the oxygen atoms in the tetragonal 
unit cell centered on site f. Values Ed and Sp are the 
copper and oxygen hole energy levels respectively, and p 
is the chemical potential. Operators and P(a)s,a de¬ 
scribe the destruction of the holes with spin a at orbitals 
dx 2 _y 2 and p(a) = {px,Py}', n’^d)f ^-nd are the cor¬ 

responding hole number operators. The hopping tp^p^, 
is between two nearest-neighbor oxygen sites, while tp^d 
corresponds to the hopping between neighboring copper 
and oxygen orbitals. The explicit form of hopping terms 
depends on the chosen phase conditior^ and is consis¬ 
tent with the paper.— The largest energy scale in the 
problem is defined by the set of the Coulomb repulsion 
parameters which consists of intra-atomic interactions Ud 
and Up and the nearest-neighbor copper-oxygen Coulomb 
parameter Vpd- 

The ab initio hopping parameters and single electron 
energies of the underlying p — d model m have been 
obtainedi^ earlier for La 2 Cu 04 using a Wannier function 
projection procedure^!. Their values are listed in TableHl 
It should be noted that parameters in this Table differ 
from the values in Table HI of Ref [l^ since we do not 
normalize them by tpd- The Coulomb parameters are 


+ E {pUg.P(c.'}s'a + ^-C-) + 

a,a',g/g',(T 

g,C7 Q:,f ,g,(7,(7' 


taken from the Ref. [1^ where the following values have 
been obtained by fitting to experimental ARPES data: 
Ud = 9, Up = 4:, and Vpd = 1.5 (all values are given in 
eV). 

The next term in Eq. ([T]) is the phonon Hamiltonian 
Hph which includes the only in-plane full-breathing oxy¬ 
gen vibrations. It is sufficient for the purpose of the pa¬ 
per although more realistic approach should take at least 
several modes with the largest coupling constant into ac¬ 
count. In cuprates these modes^Si^l involve breathing 
Cu-0 stretching motions of the planar and apical oxygen 
atoms and out-of-plane buckling motions of the planar 
oxygen atoms. Nevertheless to reasonably simplify the 
problem the phonon Hamiltonian Hph has been written 


Parameter 

£d 

£p 

tpd 

tpp 

value, eV 

0 

0.91 

1.36 

0.86 


TABLE I: Hopping parameters and single-electron energies 
of the p — d model © obtained in the framework of the 
LDA+GTB methodic and used in the present paper. 
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down as 

Hph — filOhr ^ ^ 2 ^ ' 

Operators and e\ describe the destruction and cre¬ 
ation, respectively, of the local phonon at site f with the 
frequency ujbr ■ They are linearly related to the displace¬ 
ment operator via the canonical transformation: 



that convert a quadratic form of the initial vibration 
Hamiltonian to a diagonal one. Here is the unit polar¬ 
ization vector, which is assumed to be a constant inde¬ 
pendent of reciprocal space vector q as well as the phonon 
frequency LObr ■ 

The phonon-induced renormalization of the electron 
energies is given by the Hamiltonian Hepi 

Hepi — 

f,(T I 

+ (4aPg<T + h.c^ ■ (5) 

The first diagonal term results from modulation of cop¬ 
per on-site energy by the oxygen displacements from 
their equilibrium positions. Since the value Ed increases 
with the Cu-0 bond length increasing, the index Si = 0 

for displacements and S'; = 1 for 

where and By are the unit vec¬ 
tors in the directions x and y. The second term in Eq. ® 
is off-diagonal one and describes the dependence of Cu-0 
hopping energy on the length of Cu-0 bond. 

It is convenient to define the dimensionless param¬ 
eter of electron-phonon interaction which will be used 
throughout this work as a measure of the electron-phonon 
coupling strength, 

, _ 9d{pd) 

d{pd) 2 Muj'^W' 

I 


For the chosen set of parameters the bandwidth W of 
the occupied valence band in La 2 Cu 04 is equal to 2.2 eV. 
The phonon frequency is taken as ujbr = 90 meV in accor¬ 
dance with the measured^^ value. In the present paper 
the parameters ^d{pd) ’^^.ry between 0 and 0.5 values. 

To proceed in the scheme of generalized tight-binding 
method we need to separate the total Hamiltonian ([T]) 
into the intracell and intercell parts. The concomitant 
problem of non-orthogonality of all operators related to 
the oxygen sites of adjacent cells is solved explicitly via 
the canonical transformatior^ that introduces new oper¬ 
ators in /c-space. The detailed description of this proce¬ 
dure for the Hamiltonian of the p—d model can be found 
elsewhere.— >2^ The new oxygen hole operators are a lin¬ 
ear combination of the Fourier transforms of the original 
Pix)qa and P(y)q<T orbitalsi 

^qfT ('^3:P(a:)qcr “I" ^yP{y)cia) ; 

Mq 

^qfT {^xP(a;)q(T ^yP{y)^(T^ i (^) 

Mq 

where s^^y) = sin a is a lattice parameter, and 

P<i — + Sy. In the coordinate space, the group or¬ 

bitals aqcr and are Wannier-like oxygen wave func¬ 
tions of aig and big symmetry, respectively, which are 
centered at the copper site f and spread over several 
neighboring sites. 

To diagonalize the phonon part of the Hamiltonian o 
the procedure similar to the Shastry canonical transfor¬ 
mation has been successfully applied in Ref. . Fourier 
transforms of the new phonon operators are given by 
equations 

^q “ ('5xe(a:)q T ^y^(y)q) : 

Mq 

2 

Bq = ~~ (■5ye(j,)q — Sa;e(y)q) . (8) 

Pq 

After the orthogonalization (00 the total Hamiltonian 
can be written as a sum of the individual unit cell part 
He and the intercell contribution Hec 


H = He + Hee, 


He = Y,Hfa, 

f ,(T 


He 


f.f'.lT 


Hfa =Y^^P- 9 ) + ^-c-) + 2Y^t 


/3"-03)f’^(/3)f 






p y u' 

-\- Aj. -|- -|- Aj.^ ^(d)f T ‘^\dPo (a| + A^^ T h.c.^ , 

Hff'a = — (dl^bf,^ -|- h.c.^ — 2tppi'fj, (bl^b^,^ + h.c.^ + 

+ 2A^^jj, (Af + Aj^ Tf(d)f' + “^^pdPs'h. (aJ -I- ^f) + ^-c-) ■ 


(9a) 


(9b) 

(9c) 
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Here index /3 = {d, 6} enumerates the plane orbitals of 
copper and oxygen, respectively. The terms containing 
operators of aig oxygen orbital as well as H-type phonon 
excitations are omitted in the Hamiltonian because 
of their insignificantly small influence on the low lying 
local eigenstates. Expressions for the new parameters 
and renormalizing coefficients from Eq. m are presented 
in the Appendix A. 

We should to emphasize that electron-phonon interac¬ 
tion in the derived Hamiltonian has the non-local charac¬ 
ter. The procedure of orthogonalization results in renor¬ 
malization of all matrix elements in Eq. ® which become 
strongly distance-dependent even for the initial parame¬ 
ters are taken to be non zero only for the nearest neigh¬ 
bors. 


III. THE POLARONIC VERSION OF THE GTB 
METHOD 

The GTB method has been proposed to calculate the 
band structure of compounds in the limit of strong elec¬ 
tron correlations. When Coulomb energy is much higher 
then kinetic energy the idea of bare electron as zero ap¬ 
proximation of a theory does not work. Contrary to 
the conventional tight-binding method the local states 
in GTB approach are not free electron ones but rather 
quasiparticle excitations between multielectron terms of 
d” and d"^^ configurations. 

It should be clarified what these quasiparticles are. 
The localized multielectron d” configurations for a sep¬ 
arate ion in the crystal field can be easily determined 
from local electroneutrality. Lets denote m possible for 
the given configuration d” terms as Em (n). One of them 
is the ground term Eq (n) that is occupied at zero tem¬ 
perature. An excitations from the ground to any term 
with the energy AEmo = Em (n) — Eq (n) are the local 
Bose-type quasiparticles such as excitons, magnons and 
so on. If external electron comes to the given ion the 
later changes its configuration to d"+^ with its own en¬ 
ergy spectra Em' (n -I-1). The electron addition energy 
^mm' = Em (''T' + 1) ~ Em' (n) may be considered as the 
single particle excitation between two multielectron con¬ 
figurations, with the initial state Em [n) and the final 
state Em'in + 1). The interatomic interactions trans¬ 
form these local excitation energies flmm' into the quasi¬ 
particle bands flmm' (q). 

Due to the large number of initial and final states 
there are different quasiparticles with all possible pairs 
of (m, m'). It is evident that contribution of each partic¬ 
ular |m, n) —>■ |m', n -|- 1) excitation is determined by the 
corresponding matrix element (m, n -|- 1| \m', n) of the 

electron creation operator . Moreover excitations from 
empty |m,n) to empty \m',n) states have zero spectral 
weights while their energies VLmm' are defined. Non zero 
spectral weight of the quasiparticles results from total or 
partial occupation of participating multielectron terms. 

Therefore in GTB picture the correlated electron is 


treated as a linear combination of different quasiparti¬ 
cles and more important that each of them has its own 
quasiparticle weight. It results in the crucial difference 
between free electron picture and GTB one. The follow¬ 
ing spectral weight redistribution over these quasiparti¬ 
cles defines the underlying effects of the band structure 
formation in correlated systems. 

The polaronic version of the GTB method is the 
natural development of the approach to the systems 
with strong electron-electron and electron-phonon in¬ 
teractions. While in GTB picture the quasiparticles 
are formed via the electron-electron interactions in p- 
GTB method the formation of quasiparticles is owing to 
the presence of both the electron-electron and electron- 
phonon interactions. 

IV. EXACT MULTIELECTRON AND 
MULTIPHONON EIGENSTATES OF CUO 4 
CLUSTER 

Undoped La 2 Cu 04 has a mixture of (hole on cop¬ 
per) and d}^p^ (hole on oxygen) configurations. Both are 
related with one hole per unit cell (here CUO 4 cluster). 
The electron addition results in the d}^p^ local configu¬ 
ration without holes. The electron removal results in a 
mixture of two-hole local configurations <Pp^, d^^p'^ ^ and 
That is why one has to start the GTB procedure 
with exact diagonalization of the local Hamiltonian He 
given by ([5]) in the 3 subspaces of the Hilbert space with 
the number of holes n/j = 0 , 1 , 2 . 

A. Subspace with Uh = 0 

The eigenstates can be written as 

\0,u) = \Q)\u),V = 0, I, ..., (10) 

Here |0) means the hole vacuum, corresponding to the 
electronic configuration \u) denotes muliphonon 

state with number of phonons Uph = v that results from 
V X action of phonon creation operator A~^ on vacuum 
state |0) of harmonic oscillator: 

M = ^(A+yio}. (11) 

V A 

B. Subspace with Uh = 1 

Without EPI there are two eigenstates corresponding 
to the bonding and antibonding mixtures of dx 2 -y 2 cop¬ 
per and Px,y oxygen orbitals. The i-th cluster eigenstate 
with one hole is a spin doublet with projection of spin cr, 
that may be written in the following way 

^max 

llc^u) = +<^iv\b^)W))- (12) 

1^=0 
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FIG. 1: The distribution of a hole among copper d and oxygen b orbitals as function of (a) the diagonal EPI Ad and (b) the 
phonon number v for small coupling and (d) the same for large coupling, (c) the ground state | Icr, 0) ratio of the square root of 


the average square of oxygen displacement to the lattice parameter 
dependencies on the diagonal EPI. 


^0 


polaronic shift E-po, and number of phonons {Nph)^ 


Here \dcr) = |0), \ha) = 6+ |0). In general, the ground 

and excited eigenstates m characterize the electron sur¬ 
rounded by a cloud of phonons, e.g. the polaron. The 
polaronic shift Epi of the term m may be calculated as 
a difference in the term i energy with and without EPI: 
Epi (A) = Ei (A) — Ei {X = 0)^ It is instructive also to 
calculate a square root of the average square of oxygen 
displacement 






only one summand in Eq. m with v = 0. For small 
EPI Xd = 0.01 the maximal probability is for 0-phonon 

state ^|cgo|^ + |coo|^) |0), phonon cloud around electron 
is very thin. For Ad = 0.1 the maximal probability to find 
a hole either on copper or on oxygen occurs for three- 
phonon V = i state with rather large contributions from 
1-, 2- and 4-phonon states, similar demonstration of pola¬ 
ronic effect on the local single-hole ground state has been 
obtained in the paperi^ With increasing EPI the max¬ 
imal probability shifts to the multiphonon states with 
1 / = 15 for Xd = 0.3 and = 35 for Xd = 0.5 (Fig. [TJl). 


for the cluster centered at site f in the eigenstate |nd,i). 
This value allows to restrict the EPI parameters not to 
get the lattice melting following the Lindeman criteria. 
For values Ad, Xpd <0.5 the Lindeman criteria {iE)q “C 
a is fullfilled (Fig. [!}:). 

For simplicity we start our discussion with the effect of 
diagonal EPI when the off-diagonal term Xpd = 0. The 
monotonic distribution of one hole among copper and 
oxygen orbitals as function of Ad is shown in (Fig. [TJi). 
A structure of polaronic state m is determined by 
the coefficients cf^ and that are shown in (Fig. [T]d) 
for different phonon numbers. Without EPI there is 


The maximal number of phonons iVmox in (Fig. [T^) 
is a parameter of the theory however the choice of this 
parameter depends on value of Ad(pd)- For any given 
EPI coupling we have calculated the hole distribution vs 
phonon number v for different N^ax, and increase the 
Nmax value up to stable zero contributions for higher 
phonon numbers. Thus, for Ad = 0.5 Nmax = 60, and for 
Ad = 1 Nmax = 120. For small Ad effect of EPI on the 
occupation number of d-orbital is weak because the po¬ 
laronic shift Ep is small vs covalence that is determined 
by tpd, lEpI <C tpd- Hole hopping from oxygen to copper 
takes time Ri 1/tpd and it is faster than polaron forms 





































































































6 




FIG. 2: The effect of both diagonal and off-diagonal EPI on (a) the distribution of a hole as function phonon number and Xpd, 

/ (u'^ \ 

(b) ratio of the square root of the average square of oxygen displacement to the lattice parameter * , polaronic shift Epi, 

and average number of phonons {Np|l)^, (c) the energy levels of a single hole ground and excited states, (d) the change of the 
copper and oxygen orbital occupations for the ground and excited states. 


a potential well and localized in it for the time ~ 1/Ep. 
A charge carrier is untrapped by the deformation held 
and delocalized, we can call such state a large radii local 
polaron (large local polaron). Term “local” here indi¬ 
cates the state without dispersion. The polaronic disper¬ 
sion will be considered in the next section, nevertheless 
classihcation of polarons as large and small radii corre¬ 
sponds to delocalized/localized states. The stronger EPI 
the larger is population of hole on copper orbital that 
results from the hole self-trapping, because oxygen holes 
are more mobile than copper ones. Smoothly with Xd 
growth at Xd = 0.03 — 0.04 the phonon cloud transforms 
from the narrow distribution with the maximum at 0- 
phonon component in the eigenstates (1121) to the multi¬ 
phonon components. Nevertheless the partial occupation 
of oxygen orbital takes place for large EPI also. Only in 
the non-realistic limit = 1 we have found full occupa¬ 
tion of the d-orbital (Fig. [TJi). 

Now we consider the effect of both diagonal Xd and 
off-diagonal Xpd EPI. Contrary to the diagonal EPI the 
off-diagonal one has a trend for hole delocalization. It re¬ 
sults in a partial compensation of both EPI contributions. 
For the same value of A the multiphonon contribution to 
the eigenstate (1121) is much smaller (compare Fig. [5^ and 


Fig. [T}d). From Fig. I2t i it is clear that increasing of Xpd 
up to the Xd value results in decreasing the size of the 
phonon cloud and the multiphonon weight. With param¬ 
eters XdjXpd there are too many variants, for simplicity 
we restrict ourselves further by the case Xd = Xpd = A. 
In this case we have found a critical value Ac = 0.314 
separating properties of our local polaron. At A < C Ac 
the square root of average square of displacement (u^), 
number of phonons (Ap/,), and polaronic shift Ep are less 
than for the diagonal EPI regime (compare Fig. [2 ]d and 
Fig. [TJ:). This is effect of partial compensation of diag¬ 
onal and off-diagonal EPI. At small A the redistribution 
of hole between copper and oxygen vs A is almost absent 
(Fig.[2j:), polaronic shift of the ground and excited eigen¬ 
states (|T^ is almost the same (Fig. [ 2 ) 3 ), the number of 
phonons (Nph) is close to the case of EPI absence, the 
difference (l,i -I- 1| Nph |I,i-|-l) — (l,i| Nph |I, i) is close 
to 1. These are the large local polaron states. Approach¬ 
ing the critical value Ac all characteristics of the ground 
and excited terms m have drastically changed. Non 
monotonic dependence appears for the hole occupation 
numbers (Fig. [2}:), the average square of displacement 
number of phonons {Nph)^ (Fig.j^Jo) in the ground 
and excited terms, the energies of the ground and excited 
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FIG. 3: The copper and oxygen hole distribution vs phonon 
number in the single-hole ground and three excited states for 
equal parameters of the diagonal and off-diagonal EPI (a) 
below the critical value, (b) in the critical point and (c) above 
it. 


terms are close to crossover (Fig. [2ji). 

We have compared the polaronic effect in the ground 
and several excited states in Fig. [31 Below the criti¬ 
cal value there is maximal contribution of = 0 in the 
ground state while the inultiphonon contributions domi¬ 
nate in the excited states (Fig. |3Ji). At the critical value 
A = Ac = 0.314 both 0 -phonon and multiphonon contri¬ 


butions exist in the ground state, one is typical for the 
large and the other for the small polaron(Fig. [5 }d). Si¬ 
multaneously in the first excited state (Fig. [5 }d) we have 
noticed the increase in the 0 -phonon contribution vs the 
same excited state in Fig. [3^. 

All these results indicate that at A = Ac = 0.314 
there is a crossover between large and small local polaron 
states. This critical behavior is different to the case of 
the diagonal EPI where the transformation between large 
and small local polaron states is smooth. For A = 0.32 
in the ground state with rih = 1 we have found almost 
full occupation of the copper orbital and empty oxygen 
one (Fig. |2t, Et). Further increase of A results in a small 
decreasing of copper hole occupation while the average 
square of displacement number of phonons (A^p/t)g, 

and polaronic shift |ifpo| continue to grow. 


C. Subspace with nn = 2 

Two-hole states of the Cu 04 cluster may be singlets 
or triplets. For the low energy theory only singlets are 
essential. We write down the exact eigenstates in the 
form 

N 

| 2 . 7 ) = + 

+ \did^) W) + c“ \hh) W)) (13) 

Here the first term corresponds to the Zhang-Rice (ZR) 
configuration \3dP2p^), the second one to the two holes 
on copper |3c?®2p®), and the third one has two holes on 
oxygen |3d^°2p^). The two-hole state appears in the the¬ 
ory due to electron removal process from the single-hole 
state. The additional hole occupy mainly the oxygen or¬ 
bital. From Fig. 0^ it is clear that in the absence of EPI, 
at Xd = Xpd = 0 the number of copper holes has increased 
by 0.04 in comparison to Fig.[TJi while the number of oxy¬ 
gen holes has grown almost by 1. In the ground two-hole 
state the main contribution is given by the ZR contribu¬ 
tion, the minimal weight has |d 4 ,dt) contribution, and the 
164 , 6 -t-) configuration has intermediate weight (Fig. 0 Jd). 

If electrons were free of the Coulomb interaction one 
would expect a doubling of many polaron characteris¬ 
tics vs single-hole state ones. Due to the strong electron 
correlations in cuprates situation is different. Because 
the second hole occupies mainly the oxygen orbital the 
diagonal EPI results in a small changes vs single-hole 
eigenstates. Eor Xd <0.1 the average square of displace¬ 
ment number of phonons (Nph), and polaronic 

binding energy increases less then 25% in comparison to 
the single-hole eigenstates (Pig. lit). Purther increase of 
the diagonal EPI provides a smooth growth of the cop¬ 
per hole population, at Xd = 0.24 it becomes equal to 
the oxygen hole population (Fig. |lt) and is close to 1 at 
Xd = 0.9. Thus we have found for two-hole eigenstates 
the effect of competition between the Coulomb repulsion 
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FIG. 4: The distribution of the two holes among copper d and oxygen b orbitals as a function of the diagonal EPI A and phonon 
number v for (a) small and (b) large coupling, (c) the ground state ratio of the square root of the average square of oxygen 


displacement to the lattice parameter 
diagonal EPI. 




, polaronic shift Epo, and average number of phonons (Nph)Q dependences on the 


and effective attraction mediated by EPI. Similar compe¬ 
tition has been revealed earlier by the Quantum Monte 
Carlo methodThe transformation of the phonon cloud 
with the maximum at 0-phonon component to the multi¬ 
phonon maximum for two-hole states m with diagonal 
EPI occurs smoothly in the region = 0.025 — 0.03 
(almost the same as for single-hole states dT^ '). never¬ 
theless the evolution the large local polaron to the small 
local polaron continues up to Xd = 0.9(Fig. |3Ji) 

In the regime of equal diagonal and off-diagonal EPI 
up to A < 0.314 the population of oxygen holes negligi¬ 
bly increases and the copper holes population decreases 
with strengthening the EPI(Fig. [5^). At the critical 
point A = 0.314 there is a crossover between the differ¬ 
ent two-hole configurations, the maximal population is 
acquired by the |d4.(it)-configuration. The copper holes 
number increases sharply almost by 1 (Fig. [S^) at the 
critical point. Corresponding jumps have been revealed 
for the average square of displacement and num¬ 

ber of phonons Nmax (Fig. Eb). The maximal number 
of phonons sharply changes from = 15 till r/ = 190 
(Fig.Et). 

Summarizing this chapter we plot schematically in 


Fig. inii the general structure of the relevant Hilbert 
space with three subspaces with the number of holes 
rifi = 0,1,2. With maximal phonon number each 

electronic level is splitted into (A^max + 1) sublevels. In 
the hole vacuum subspace there are (A^max + 1) equidis¬ 
tant sublevels numerated by number of phonons and 
separated by a phonon energy hubr- In the single-hole 
subspace there are two blocks of spin doublets with the 
2 (A^max + 1) sublevels each, corresponding to the bond¬ 
ing and antibonding hole orbitals. In the two-hole sub¬ 
space there are three singlet blocks of (A^max + 1) sub- 
levels and the triplet block of 3 (A^max + !)■ In subspaces 
Uh = 1 and Uh = 2 the energy levels are not equidistant 
and the phonon number is not a quantum number. 

Due to the long-range antiferromagnetic order in the 
undoped La2Cu04 with two sublattices A and B the lo¬ 
cal hole states from Eq. (HID are subject of the effective 

.2 

exchange interaction J ~ that results in the split¬ 
ting of states |Icr, i) and |1 ct,*), here a means —a. From 
the chemical potential equation for undoped La2Cu04 
at T = 0 we found the only occupied term, that is the 
= 1 ground state marked by cross in Fig. Ac¬ 
cording to the GTB approach (see Section Hill above 1 the 
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form the valence (v) band. Similarly, the electron addi¬ 
tion to the occupied single-hole term results in formation 
of a set of hole vacuum multiphonon states (fTOll and the 
conductivity (c) band. 


V. POLARONIC BAND STRUCTURE IN THE 
GENERALIZED TIGHT BINDING METHOD 

The mathematical tool that allows to work with 
multilevel lattice system with a set of orthogonal 
and normalized local eigenstates {|p)} = {|n.i)} = 
{|0,i^), |lcr,i), |2, j)} is given by the Hubbard operators 
X'P'^ = \p) (glf^ Each diagonal operator X'pp determines 
the occupation of the eigenstate |p), while the nondiag¬ 
onal operator Xf^ describes the excitation at the site f 
from the initial state jg) to the final state \p). If the 
change of electric charge during the excitation is odd 
this excitation is the Fermi-type quasiparticle. According 
to the definition of our inultielectron and inultiphonon 
eigenstates, the Hubbard fermion in p-GTB is a polaron 
and we call it the Hubbard polaron. Due to completeness 
of the local set of eigenstates {|p)} each single-hole anni¬ 
hilation operator at site f is given exactly by the linear 
combination of the Hubbard fermions 

We can write down the hole annihilation operators on 
corresponding orbitals /3 = 6, d as the linear combinations 
of the Hubbard fermions 



0 100 200 300 400 

Phonon number v 


(14) 

pq 

The matrix elements 7 o.(/ 3 ) (?"?) = {p\ k) for orbital 
/3 are calculated straightforwardly because we know all 
eigenstates \p) and \q). Quasiparticle transitions between 
states within one Hilbert space sector are described by 
Bose-type Hubbard operators to distingwish Bose 

operator from Fermi one we use for the former the nota¬ 
tion Z. 

The phonon number is not a quantum number in the 
single-hole m and two-hole m states therefore the 
phonon annihilation operator on the site f is given by 


FIG. 5: The effect of the equal diagonal and off-diagonal EPI 
on the distribution of the two holes as the function of (a) 
EPI and (c) phonon number, (b) the ground state ratio of the 
square root of the average square of oxygen displacement to 


the lattice parameter ^ ^ ° , polaronic shift Epo , and average 
number {Nph)Q of phonons. 


electron removal from this initial state results in forma¬ 
tion of the set of hole quasiparticles with different final 
states in the subspace Uh = 2. These excitations are 
shown schematically by arrow in Fig. [5^. Their disper¬ 
sion is studied in the next Section |Vl all these excitations 


Af = (15) 

n—0 ij 

where ja [ni, nj) = (m| Af \nj). 

Since we know all coefficients in Eq. m and (USD it 
is easily to write down the total Hamiltonian He -h Hec 
from Eq. [5] in terms of Hubbard operators, then it con¬ 
sists of two parts. The first one, Hei-int = He + Hfe~^\ 
is similar to the Hubbard model but constructed on oper¬ 
ators describing Hubbard polaron excitations instead of 
usual Hubbard fermions ones. The second part, Hei-ph = 
Hfe ~^^, results from EPI in initial Hamiltonian and con¬ 
tains Bose excitations. 
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FIG. 6: (a) The schematic picture of the multielectron and multiphonon local eigenstates for number of holes nu = 0,1, 2 
per Cu 04 unit cell and maximal number of phonons Nmax = 2. In the two-hole sector, the irrelevant triplet states with 
higher energy are not shown. Arrows between subspaces (0,1) and (1,2) schematically indicate polaronic quasiparticles, (b) 
Quasiparticles dispersion without EPI. Solid lines correspond to the conductivity and the valence bands of electrons in the 
antiferromagnetic phase. Dotted horizontal lines correspond to the dispersionless Franck-Condon resonances with zero spectral 
weight in the absence of EPI. 


= Hr 


TJ U _ Tjel — el , Tjel—ph 

tlcc^ Flee — Flee + ^cc ’ 


(16a) 


f n—O i 


(16b) 


= - Y. E { {llip) (0*, l^j) + lY) (isfc, 21) X 

ii' jj' kk'll' 

(lcr{b) (0*', Icrj') X^J- -b 7 ^(b) (Id-fc', 21') ^ -b /i.c.j , (16c) 

= E E E (m, nj) (m, nj) x 

f.f',(7 n,i,j k^l,r 

(|7f7(d) + |7a(d) (lCTfc,2r)|^Xf,''’^'') + 

+ E E E 2Apdp4'h (7a(A) +7<t(A) X 

crn,z,_7 k,k',l,l' 

{ (7:(d) (Ok, lak') + 7 :(,) (la;, 21') x 

( 7 .( 6 ) (o;, la;') + 7 ,(,) (la^, 2k') + h.c.} . (16d) 


Here the intracluster term He contains only energies Eni 
of the local eigenstates of the cluster with diagonal op¬ 
erators X™’”*, where n = 0,1,2, and index i enumer¬ 
ates the polaronic eigenstate in the Hilbert space sector 
with hole number n. The intercluster terms result in 
the polaronic hopping and dispersion. Due to the inter¬ 
cluster contribution from the diagonal and off-diagonal 
EPI the polaron scattering on the bosonic exci¬ 

tations occurs. In the intercluster matrix elements the 


single electron hopping t^d, tpp and EPI A^, \pd parame¬ 
ters are strongly suppressed by the matrix elements from 
Eq. (imi and (fTSl) and structural factors /ifg, J^fg, Pfgh- 
These factors are strongly decreasing with distance, and 
usually only contributions from the first three neighbors 
are enough, that is why the t —f' —t"-tight binding fitting 
is rather successful in cuprates 

To obtain dispersion of quasiparticle excitations we use 
the method of equation of motion for the two-sublattice 
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FIG. 7: Effect of the diagonal EPI on the polaronic band structure for undoped antiferromagnetic La2Cu04 a,tT=lQK for 
different values of EPI parameters: (a) Ad = 0.1, (b) Ad = 0.2, (c) Ad = 0.3, (d) Ad = 0.4. The line intensity is proportional 
to the quasiparticle spectral weight. We emphasize the decreasing of the intensity scale for the conductivity band with EPI 
growth. 



FIG. 8: The effective mass of the first removal state at the 
top of the valence band for different values of diagonal EPI, 
Apd ~ 0. 


matrix Green function of polarons D'^9 {uv,v'u') = 
(^{xfQ Due to the large number of fermionic 

quasiparticles it is convenient to introduce the matrix 
Green function Dffi = {D^^}, where row m and col¬ 


umn n indexes numerate the different quasiparticles 
{p,q)n,{m) ^ n{m). Indeed, the number of quasiparti¬ 
cles (p, q) is finite, so each may be enumerated just by 
the number n {m) which has a meaning of the quasipar¬ 
ticle band index. Total Green function in the matrix form 
looks like 


Dff' 


D§P D§P ) 


(17) 


For the Hamiltonian with two-particle electron-electron 
interactions the generalized Dyson equation for the ma¬ 
trix Green function^^ in the Fourier transformation reads 


Dk (w) 


^0 ^ (^) ~ (^) + Sk (to) 


Ac(w) ( 18 ) 


Here G'(j'^(a;) is a local propagator determined by the 
multielectron eigenstates |p) and |q), = 'y{m)^{n)t]^ 

is the intersite hopping matrix, and <k is a bare disper¬ 
sion. Besides the self-energy E(k, w) the unusual strength 
operator P(k, w) appears in Eq. (fTSl) . It results in the 
redistribution of the QP spectral weight (or oscillator 
strength) and in renormalization of quasiparticle disper¬ 
sion, that are the intrinsic features of SCES. Recently the 
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k 

FIG. 9: Effect of equal diagonal and off-diagonal EPI on the polaron band structure of the undoped antiferromagnetic La2Cu04 
at T — IQ K for different values (a) A = 0.1, (b) A = 0.2, (c) A = 0.3, (d) A = 0.35. 



FIG. 10: The effective mass at the top of the valence band at 
(f ’ f) fiffual diagonal and off-diagonal EPI. 


ARPES line shape had been discussed with Eq. (IT^ and 
the odd in (k — lipermi) contribution to the momentum 
distribution curve had been found due to the imaginary 
part of the strength operatori^ 

The simplest nontrivial solution (Hubbard I approxi¬ 
mation) is usually used in the cluster perturbation the- 
or y34,35 ^ It can be obtained from Eq. (1151) when Ek (w) = 


0 and Pk{uj) = 5pqF{pq). The so-called filling factor 
F {pq) is given by the sum of the initial and final state oc¬ 
cupation numbers Ff (pq) = {Xf^) + {Xf'^) and is strongly 
dependent on sublattice magnetization, doping, and tem¬ 
perature. The appearance of the filling factor is very im¬ 
portant difference of the Hubbard fermions vs bare elec¬ 
trons. In particular, this factor makes irrelevant many 
excitations from empty \p) to empty jg) states with de¬ 
termined energy Ep — Eq. 

In the Hubbard I approximation we have obtained the 
following Dyson equation for the polaronic matrix Green 
function 


(19) 

Here = F {pq)/{uj — D {pq)), D {pq) is the local quasi¬ 
particle energy of the multielectron and multiphonon 
eigenstates and is the matrix of combined p — d and 
p — p hopping that provides the band dispersion Wk in 
the absence of EPIi^^ 

Matrix contains terms of intercluster EPI 

^{pd)k ^(pd)k ""^hich are presented in the Ap¬ 
pendix m- For the typical EPI parameters we have 
estimated the ratio ^ 0.01, and the terms 
and AI^p^)k Provide small contribution to the dispersion. 
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Thus, the local intracluster effects of EPI discussed in 
the Section IIVI are the most important for the polaronic 
bands formation. In the absence of EPI the only po¬ 
laronic quasiparticles with non-zero spectral weight are 
possible between the multielectron n, n -I- I terms with 
equal number of phonons. At zero temperature this con¬ 
dition results only in 0 — 0 Franck-Condon resonances 
between terms with the phonon number u = Q. These 
Hubbard fermions dispesion is shown in Fig.|6}D by a solid 
line with different spectral weight in various parts of the 
Brilloin zone due to long range antiferromagnetic order. 
All other polaronic excitations corresponds to the dis¬ 
persionless Franck-Condon resonances with zero spectral 
weight in the absence of EPI (Fig. [61 d). 


In general, the EPI results in the hybridization of the 
Hubbard subbands and the Franck-Condon resonances. 
The low energy part of the polaronic bands is shown in 
Fig. [7] and Fig. O Effect of the diagonal and both diag¬ 
onal and off-diagonal EPI on the polaronic band struc¬ 
ture is different. Weak diagonal EPI introduced in ([9]) 
mainly modifies the conductivity band while equal di¬ 
agonal and off-diagonal coupling more strongly effect on 
the valence band. Thus at A = 0.1 we found the hy¬ 
bridization splitting and finite spectral weight for sev¬ 
eral Franck-Condon resonances in the conductivity band 
(Fig.[^). The number of splitted subbands and the value 
of the minigaps between these subbands becomes larger 
for A = 0.2 (Fig.jTb). For the valence band there are also 
splitted subbands separated by smaller minigaps. The 
top of the valence band, particularly the first removal 
state at the (7r/2,7r/2), is unaffected by small diagonal 
EPI. For larger EPI the spectral weight of the conduc¬ 
tivity band is strongly suppressed, it is transferred to the 
higher energy bands shown in Fig. [6^. For \d = 0.4 the 
first removal state is also strongly renormalized (Fig.[7jl). 
The effective mass of the hole in the first removal state 
sharply increases above Xd = 0.35 while for Xd < 0.35 
the effective mass my along the (0,0 ) — (tt, tt) direction 
and m*^ along the (tt, 0 ) — (0, tt) direction increase rather 
weakly with EPI growth (Fig. [8]). This sharp increase of 
the effective mass results from the polaron autolocaliza¬ 
tion, similar effect has been obtained in Refs. [^IT§|. 


In the regime Xd = Xpd = X both the conductiv¬ 
ity band and the middle part of the valence band are 
strongly renormalized by the subbands splitting and 
spectral weight redistribution (Fig.|9K-d). In this regime 
we have found the critical coupling value changing in the 
properties of local single- and two-hole states (see Fig. [3] 
and Fig. [5] above). The sharp increases of the effective 
mass occurs also at A = Ac iFig. ITUl) . that is the manifes¬ 
tation of the polaron autolocalization. 


VI. POLARONIC SPECTRAL FUNCTION 

The polaronic spectral function is given by 
^(k,w) (k,w) = 

E 7a(/3)(m)7*(/3)Hlni^k’”(‘^ + *^)- 

(T(3mn 

Splitting of the Hubbard fermion bands to the many hy¬ 
bridized subbands results in the series of narrow peaks 
for a given wave number. Each peak results from some 
Franck-Condon resonance and is related to the multi¬ 
phonon excitation. For Xd < 0.2 the main peak of the 
first removal state at the top of the valence band with 
k = (f, f ) corresponds to the 0 — 0 resonance between 
the 0-phonon single- and two-hole ground states from 
Fig. |6^. The multiphonon contributions to the spectral 
function are negligibly weak and shifted down in energy 
(Fig. [TThl. With the diagonal EPI increasing the 0 — 0 
peak decreases while the set of multiphonon peaks ap¬ 
pears fFig. [TTbl. For the eqnal diagonal and off-diagonal 
EPI situation is qualitatively similar (Fig. [TTh.db The 
shift of the multiphonon peak from the 0-0 resonance 
with almost zero spectral weight is larger in Fig. [TTH . 

The effect of the quasiparticle finite lifetime on the 
spectral function is also shown in Fig. Illb .d. It is mod¬ 
elling by the different Lorenzian width 5. With increas¬ 
ing 5 we reproduce the formation of one wide peak in the 
spectral fnnction from the sum of several Franck-Condon 
resonances. This mechanism of the large linewidth in the 
ARPES experiments in the nndoped cuprate has been 
discussed in the paperi^ The large shift of the spectral 
intensity below the nominal top of the valence band in 
the absence of EPI given by 0 — 0 resonance has been also 
found in ARPES measurements)^ 


VII. DISCUSSION OF THE RESULTS 

In this paper we have developed the general approach 
to the electronic strncture of Mott-Hubbard insulators 
with strong electron correlations and strong electron- 
phonon coupling. The polaronic version of GTB (or p- 
GTB) method is a variant of the cluster perturbation 
theory with exact diagonalization of the intracell part of 
the Hamiltonian, redefinition of all local fermionic and 
bosonic operators as the linear combination of the quasi¬ 
particles that are excitations between mnltielectron and 
multiphonon initial and final states. The Fermi type 
excitations are Hubbard polarons. The formally exact 
generalized Dyson equation has been solved in the con¬ 
ventional for the cluster perturbation theory Hubbard I 
approximation that results in the polaronic band struc- 
tnre of hybridized Hubbard fermions and local multi¬ 
phonon Franck-Kondon resonances. The polaronic spec¬ 
tral weight is strongly dependent on the value of the EPI. 
We have carried out all calculations here within the three 
band p — d model. It is evident that the p-GTB may be 
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CO, eV CO, eV 


FIG. 11: The spectral function of the first removal state at for Xd = 0.2 (a), and Ad = 0.4 with two values of line 

width 5 — 0.001 eV and 5 = 0.05 eV (b). The same for Ad = Xpd ~ A with A = 0.2 (c) and A = 0.35 (d). The multiphonon 
Franck-Condon resonances are almost absent for delocalized polaron at (a),(c) and appears for localized polaron at (b),(d). 


straightforwardly generalized to the multiband realistic 
model with all copper d-orbitals and all oxygen p-orbitals, 
and to the more realistic treatment of the phonon sys¬ 
tem. The restriction of the Hubbard I approximation 
is not obligatory for the GTB method, for example the 
self-consistent calculation of the self-energy and spin cor¬ 
relation functions for the t — J model had been carried 
out in the X-operators technique in the spin liquid state 
of the doped cuprates^ and within non-crossing approx¬ 
imation!^ 

The band dispersion and the spectral weight of Hub¬ 
bard polarons are strongly temperature and doping de¬ 
pendent. The finite temperature generalization of the 
GTB band structure calculations has been demonstrated 
recently for the undoped La 2 Cu 04 where the insulator 
state is shown to exist both in the antiferromagnetic 
phase below the Neel temperature and in the paramag¬ 
netic phase above the Neel temperature!^ Due to the 
finite volume of this paper we have restricted ourselves 
here to only undoped cuprate, the discussion of the ef¬ 
fects of doping and finite temperature on the polaronic 
band structure will be given in a forthcoming paper. 

We have shown that in general polaron is character¬ 
ized by a broad distribution of the phonon numbers in 
surrounding cloud. Depending on the EPI coupling the 


maxima in the phonon distribution are given by the 0- 
phonon or multiphonon contributions. The former corre¬ 
sponds to the delocalized large polaron while the latter 
describes the localized small polaron with the crossover 
from large to small polaron with increasing EPI. Previ¬ 
ously similar behavior has been found by the diagram¬ 
matic Monte Carlo method^ within the t — J-Holstein 
model. The new results in our p-GTB approach have 
been obtained when both diagonal and off-diagonal EPI 
were considered. We have found their partial compen¬ 
sation in formation the multiphonon cloud when both 
coupling parameters were equal and the sharp transition 
from the large to small polaron at the critical value of 
EPI coupling. This transition is accompanied by the po¬ 
laron localization and its effective mass divergence. The 
other our new result in comparison to the t — J Holstein 
model treatment by preceding authors is the simultane¬ 
ous transformation of both the valence and conductivity 
bands under EPI growth and different effect of the diag¬ 
onal and off-diagonal EPI on these bands. Eor the un¬ 
doped cuprate we have reproduced the large width of the 
ARPES line measured for the first removal state related 
with the loss of the spectral weight by the 0 — 0 Eranck- 
Kondon resonance and shift the spectral weight max¬ 
imum down in energy to the multiphonon resonances. 
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these features have been found experimentally in the un¬ 
doped Sr 2 Cu 02 Cl 2 J^‘^ 
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“ P^-g,g-h“ 1/A^^ X) 1/MkMqX 

kq 

[sin (A:a;/2) sin (fe/2) cos ((/ca; + fe)/2) -k (A6) 
sin {ky/2) sin (%/2) cos {{ky + %)/2)] x 

g-ik(f-g)g-iq(g-h) ^ 

Their values in depend on distance are given in Table HH 


Appendix A: Renormalized parameters of the initial 
Hamiltonian and structural factors 


The link between initial parameters and new matrix 


elements in Eq. ([S]) are given by: 


£b - £p 2tpph'QQ J 

(Al) 

Ub = bpToOOO, 

(A2) 

Vpd = lpd$ooo- 

(A3) 


Coefficients r'oo, 'I'oooo, ‘hooo are values of structural 


Appendix B: Matrices of intercluster 
electron-phonon interaction 

Matrix of intercluster EPI in the Dyson equation (fT^ 
includes two terms: 

^k'^" = <])k + <ik. (Bl) 

where -^(pd)k “ niatrices of off-diagonal interclus¬ 

ter EPI. 


factors r^fg, $ijk for single cluster, i/qq = 0.727, 

'I'oooo = 0.2109, $000 = 0.918. Values of 'kyid, $ijk are 
strongly decreased with distance increasing. 

The structural factors /ifg, ^'fg, are defined by re¬ 
lations: 


•^pdk = 25p(iCdPikX 

E 1a (pp) iZ lb 

P 

with coefficient pik = E 

gh 


(B2) 


Mfg = l/iV^/ike-'‘^(f-®), (A 4 ) 

(o') 

^pd k = 25p(iCdP2kX 

E 1a iPP) ih lb 

t'fg = l/N^{2sm{ka:/2)sm{ky/2)/pkfe~'^^‘-^~^\ ^ 

k 

(A 5 ) with coefficient p 2 k = E Po 
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